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Abstract 

In this paper we will study the existence of fundamental solutions for 

the explicit and implicit backward time dependent Schodinger equation, 

via discrete Fourier transform and its symbol for the Laplace operator. In 

both cases we will prove that the discrete fundamental solutions obtained 

►^ , converges to the continuous fundamental solution in the Zi— norm sense. 
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(3 ; 1 Introduction 

The potential theory is a useful tool to solve boundary value problems by the 
^^ , help of integral equations on the boundary. For constructive analytical consid- 

^ ' erations and also in the case of numerical applications it is necessary to have an 

explicit expression of the fundamental solution. However, the integral represen- 
tations obtained using potential theory are not suitable for an explicit compu- 
tation of the solutions, due to unacceptable convergence rates of the integral's 
numerical approximations (for more details see j7]). 

A more acceptable alternative in computational terms is given by the use of 
finite difference approximations. In fact, the connection between the potential 
theory and finite difference theory is based in the possibility of obtain explicit 
expressions for the fundamental solutions. This second approach has been ap- 
plied to elliptic operators, for example, in [3] and [8], where elliptic difference 
operators were studied and representations for their fundamental solutions were 
given. However, this approach is restricted to the stationary case. In fact, when 



one takes the time evolution into account, for example, when the operators are 
parabolic, additional difficulties appear. One of the main difficulties lyes in its 
fundamental solutions, which must be study in distributional sense. 

This paper is based on the work developed in [7], where an orthogonal de- 
composition of the function spaces in terms of subspaces of null-solutions of the 
correspondent Dirac operator were given. In order to extend this approach to 
the nonstationary case we introduce Witt basis, as done in [1]. This will gives 
us the possibility of apply techniques of the elliptic function theory to parabolic 
domains. 

For the study of the continuous Schrodinger equation 

i-A,-idt)u{x,t) = S{x,t), (1) 

we need to study the fundamental solution of the difference equation 

{{-Ah ~idr) Eh,T){hm,Tk) = 5h,T{hrn,Tk). 

For that purpose, two approaches are available to us: either one considers 
the explicit equation 

r ih if (/^™,Tfc) = (0,0) 

{{-Ah~ldr)Eh,r){hm,Tk) = I , (2) 

[ if {hm, Tk) ^ (0, 0) 

where A^, dr and 5h,T are the discrete operators of the Laplacian, the discrete 
time derivation and the discrete delta function; or we consider the implicit 
version of the equation ([2]), with an additional time-step. 

The paper will be divided as follows: on Section 2, we present some basic 
notions about Clifford analysis and finite diferences approximations. 

On Section 3 we shall study the existence and behavior of fundamental so- 
lutions for the explicit non-stationary Schrodinger equation ([2|); this study will 
be divided in three steps: first, we shall prove convergence of the discrete fun- 
damental solution to the continuous one on a bounded domain. On the second 
step, we prove convergence on a domain of type Gh x [Tq, -|-oo[, Tq > 0. We will 
finalize this section with the proof of convergence for the domain R^x]0, +oo[t-. 

On Section 4 we prove the existence of a discrete fundamental solution for 
the implicit Schrodinger equation and that this fundamental solution converges 
to the continuous one. 

In Sections 3 and 4 we will omit a general discussion in discrete spaces of 
distributions because, on the one hand we want to underline the analogy with 
the continuous case. On the other hand we are interested in convergence results 
in norms as strong as possible. Therefore, we investigate if E^.t belongs to the 
space ;i(R^ xR+). 

We will present an explicit expression for discrete fundamental solution of 
both the explicit and implicit backward Schrodinger equation and we will prove 
that these solutions belongs to the space 4°^(IK.h x Rr )■ Also, we will prove that 
these fundamental solutions converges to the continuous one in Zi— norm, when 
the size of the mesh tends to zero. 



2 Preliminaries 

2.1 Clifford Analysis 

Consider the n-dimensional vector space M" endowed with a standard orthonor- 
mal basis {ei,--- ,e„} and satisfying the multiphcation rules CiCj + e^Ci = 

We define the universal Clifford algebra C£o,n as the 2"-diniensional associa- 
tive algebra with basis given by 60 — f , as the scalar unit, and e^ = e^^ ■ ■ ■ Ch^., 
where A = {hi,. . .,/ife} C A^ = {1, . . . ,n}, for 1 < /ii < • • • < /ife < n. Each 
element x G C^o,n will be represented by x = J2a ^a^a and each non-zero 
vector has their multiplicative inverse given by j^. We denote the (Clifford) 

conjugation x — > x "" by means of its action on the basis elements 

We introduce the complexified Clifford algebra Cin as the tensorial product 
C (g) C£o,n = ^ ^ = X! ^asa, ZAeC,Aci nI , 

where the imaginary unit interact with the basis elements as iej = Cji, j = 
l,...,n. 

Finally we introduce the conjugate oi w = J2a ^aba b,sW — ^^ 'zA^'eA"^'''"' ■ 

Let now il C M" x R+ denote a bounded domain with sufficiently smooth 
boundary F = cTJ. A function u : fi — )■ C¥„ has a representation u = ^^ uasa 
with C— valued components ua- Properties such as continuity will be understood 
component- wisely. In the following we will use the short notation Lp{il), C''{fi), 
etc., instead of Lp{U, C£„), C''{i^, C£n)- For more details see [3]. 

We consider the Dirac operator D = J2^=i ^j a§" which has the property of 
factorizing the n-dimensional Laplacian, that is, D^u = —Au. A C^„-valued 
function u defined on an open domain G, is said to be left-monogenic if it satisfies 
Dm = on G. 

Taking into account [1], we will imbed R" into IR"+-^. For that purpose we 
add two new basis elements f and f^ satisfying 

f^ft'^0, fft + ftf=l, fe,+e,f=fte,+e,ft=0,j = l,---,n. (3) 

The set {f, f^} is said to be a Witt basis for R'^ and it will allows us to create a 
suitable factorization of the Schrodinger operator where only partial derivatives 
are used. 

2.2 Finite Difference Calculus 

For each 7^ G C R", (n > I), we denote the associated discrete domain Gh, 
for a fixed mesh size /i > 0, as 

Gh = {hm = (hrrii, ..., hrun) € G : m G Z"} . 



For a continuous function / on G, we denote by Rhf its restriction to the 
lattice Gh of mesh size h > 0. 

In the following we consider functions restricted to the lattices IR|, M+. We 
define the corresponding discrete ^i-spaces in the usual way: 



ueh{Rl) ^ \\u\\i,imi) 



E 



\u{hrn)\h^ < oo 



vehiRt) ^ ll«ll;i(R+)=El^(^^)l^ <°^- 



fcSN 



For a discrete function u : 



C* 



C (Xi H, given as u{hrn, kr) 



{u^ , u^ , u^ , u^) , we have the finite difference approximation for the stationary 
Dirac operators given by 



/ -9-.^ 



d;:^^ = 



■'h 

%■ 



C- -9.~!-M 



dlu' + diy 



( -d'. 



dl'u^ 



uD: 



->0 '3 



\ d-^u'^ + dW-dW J 



where 



al .0 a-3.2 , a~2..3 



DV^ = 



uDt 



dLu" - d, 



+ d: 



d^^w> 



(u^ (hm ± hcs, kr) — u^ (hm, kr)) 



dlu° - d-\^ + d-\^ 



, j=0,l,2,3, s = l,2,3, 



represent the spatial forward/backward difference operators. We remark that 
these difference operators factorize the discrete Laplacian, in the sense that 



DT-DT^u =D7^D 



+ n+- 



= -AjiU Ia = 



:ar5,>M/4, 



where Zj is the 4x4 identity matrix. In what follows we denote Df^^ and Dy^ 
as forward/backward difference Dirac operators. 

We also have the following (forward) time difference operator (for more de- 
tails see [7] and [5]) 



drU-' {hm, kr) 



u^ {hm, T{k + 1)) — u^ [hm, rk) 



, j=0,---,3. 



Further, we use the notations 
for the discrete Delta function. 



7 if fc = 
if ky^O 



We will consider the discrete Fourier transform introduced by Stummel (see 
[7] for more details) with respect to x, 

yj u{hrn, •) exp{ihrn£) for ^ £ Qh 



otherwise 

where Q„ = {c - (^1,6,6) e «' : - J < 6,6,6 < +J} ■ 

3 Explicit difference equation 

As it was indicated before, in this section we study the explicit equation (p|). 
This will be done in two steps. In the first one we study the existence of a discrete 
symbol of the operator in ^, which we use to construct Eh^r, a fundamental 
solution for discrete explicit Schrodinger backward time-dependent operator. In 
the second part we will estimate the norm 

which will allows to calculate its limit when h and t tend to zero. Here, E is a 
fundamental solution for the continuous Schrodinger backward time-dependent 
operator, where RrRhE its restriction to the lattice. 

3.1 Discrete Symbol of the Fundamental Solution 

Let us consider the equation ^. In order to simplify the resolution of this 
equation, we need to introduce the abbreviation for the symbol of the discrete 
Laplace operator 

Applying the discrete Fourier transform to ([2]), we get the equation 

{{d^ThEh,r-^^rTnE^,^r)){Lt) = ^—5^{t)^,^{^^)^ 

(27rj2 

with Xr being the characteristic function of Q/i, which has the solution 

{THEH,r){^,t) = ±.H{t){l~tTd'y-\^iO. (4) 



Using the restriction of the continuous inverse Fourier transform J^, to the 



J-f — Rh^, which acts as an inverse for J^h, we obtain 



Eh.r{hm,Tk) = iH{Tk)[{l-iT^hf ^Shjihm), (5) 

which is a fundamental solution of the discrete explicit Schrodinger backward 
time-dependent equation. 



3.2 Convergence Result 

As it was indicated before, we will now estimate the norm 

This estimation will be done in two parts. In the first we will consider the 
case of the time limited interval, i.e., we will estimate our norm in G/iX]0, TqJt-, 
with To e R"*". In the second part we will consider our norm in the G/iXjTo; +00]. 

3.2.1 Case of the Limited Time Interval 

Initially we shall study now the behavior of 

\\Eh.r ~ i?T-R/i£^||/i(GhX]0,To]x)- 

For this purpose we rewrite the equation Q in the form 

Eh,rQim,T{k + l)) 

1 + -p- 1 Eh,T{hm,Tk) 

IT 

+ —- [Eh.rihimi + 1), hm2, hm^, rk) + Ehr{h{mi - 1), ft,m2, hm3,Tk) 

+ EhT-{hmi, h{m2 + 1), hm^i^ rk) + Eh^r{hmi, h{m2 — 1), hm^, rk) 
+ Eh,T{hmi,h7n2, h{m-i + l),Tfc) + Eii,r{hmi, hm2,h{m3 — l),Tk)] , 
which implies that -E^.r is supported in a cone. From ^ and ^ we get 

\\EhA-,r)\\i,^Ml) < 1 \\Eh,r{;Tk)\\i^^j,3j < 1 p„^,(., r(fc + 1))||,^(k3) < 1. 

Further, let Tq — ttoq, with toq G N and To e IR+. By addition with respect 
to i, we get the estimation 

\\Eh,T\\h{GhX]0,TQ]r) < ll^/i,r||;i(R3x]0,To]0 - z2^ ^ '^°- (^) 

fc=l 

Now we consider the continuous fundamental solution (see [2] for more de- 
tails) 

of the continuous Schrodinger backward time-dependent operator. 



We get 

\\RhE{;t)\\,^^G,) 



(4i7ri)3 



< Hit) J2 

hmeGh 



exp 



4i 



(47rt)i 



^o/(G,.), 



where yo/(G'/i) = y. ^^ ■ Furthermore 

117? 7? FN - V ^»^(G^) - 



:^ (47rTfc)^ 



(47rT)i T 



mo -. 



fc=l 



Vo/(G,,) 



To 

(47r)3 



if 



T> 1 



< 



From (m and JHl) we conclude 



ll-E/j,^ - RrRhE\\i-^i^GHX]0,Ta]^) 



Vol{Gh) '^^^ if < T < 1 

(47r)2 



To + yo/(G,,) -^ if r > 1 

(47r)2 



(8) 



< 



To + Vol{Gh) 



^^i±> if < r < 1 

(471-)^ 



(9.) 



The previous inequality describes the approximation error of the fundamen- 
tal solution ([5]) for small values of time variable. 



3.2.2 Case of the Unlimited Time Interval 

In the following, for t — rk E M+ with t > Tq, we study 

In order to guarantee the convergence of some series and integrals, we need 
to consider the following regularized fundamental solution of the Schrodinger 
operator 

h, [x, t) = 3- exp 



(4i7ri)t 



At 



which converges, in Lp{G x K"*"), with 1 < p < +00, to the continuous funda- 
mental solution ([7]). 
We have 

\\EhA;rk) - RrRhEi-,Tk)\\i^^G,) 

= Y^ \Eh^r{hrn,Tk)- RrRhE{hm,Tk)\ h^ 

hrn^Gh 

< Vol{Gh) max \Eh_r{hrn^Tk) — RrRhE{hrn,rk)\ 

hrn^Gh 

= Vol{Gh) max \{RhTThEh,r)(hm,Tk) - {RrRhTT-^E)(hm,Tk)\ 

hrnGGfi 

= Vol{Gh) max \URhTThEhr)(hm,Tk) ~ {RrRhTT^^E''){hm,Tk)\ 

hrn^Gh 

+ \{RrRh^J'^^E^){hm,Tk) - {RrRhTT^'^ E){hm,Tk)\] 
^ I [{ThEh.r){x, Tk) - {RrT-^E){x, Tk)] exp {-ixO d^ 



< Vol{Gh) 



< 



VoljGh) 
27r 



+ max \{RrRh^^ E^){hTn^Tk)-{RrRhTT^^E){hrri,Tk)\ 

hm^Gh. 

[ \{ThEh,r){x,Tk)-{RrT-^E){x,Tk)\ |exp(-zxO|d^ 



(10) 



max URrRhJ'T-^ E^){hm,Tk) - {RrRhJ'^^^E){hrn,Tk)\ 

hmeGh 



27r 



\\{ThEh,r){;Tk) - iRrT-'E')i-,Tk)\\L,iM^) 



ii) 



■ max \{RrRh:FT-'^E''){hTn,Tk) - {RrRhTT^^E){hm,Tk)\ 

hrn^Gh 



{II) 

By the convergence of the regularized fundamental solution E"^ to the con- 
tinuous one E, we conclude immediately that the term (//) converges to zero, 
as e goes to 0+ . This fact implies that the study of pO)) depends on the analysis 
of the term (/) . In this case, we need to split this study into the outside of Qh 
and its inside. 



In R^ \ Qh, we have for the term (/) 
\\{ThEh.r)i-,rk) - (i?,^-iS^)(.,Tfc)|U,(R3) 



Hirk) {2Tkf 
(47rt) " 



2 e 2 



exp 



rk\^\ 



-) exp (-. (rfcieP - J 



iffT/c) 2 el / rfcTT^ 

exp 



In Q^ we have for (/) the estimation 

H{Tk) 



{2^)\ 



l-.rrfY-f^V 



\ e I 



exp 



rfclCP 



(11) 



ii(R3\Qh) 



exp ^« T 



-fcl^P - ^ 



ii(R3) 



g(Tfc) /2TfcY 



(^) ' (1 - -^^)'= - exp {-'-^] exp f-. frfclel^ - ^ 



Li(R3) 



g(Tfc) /2TfcY 

" (2^)i v~y 



(2^)'<'-"''^ 



'2\*; 



exp 



exp ( — z ( rkd — 



exp 



(777) 



I exp — i Tfcd 



' exp 



— 1 exp(-i(Tfc|C| - — 



L !(»■■') 



(12) 



(/y) 



From 


















e 


= 




+ sin^ 


V 2 J 


+ sin^ 


m)i 


1 

e 




< 
<- - 




1 

e 













12e 



we get 



exp,-^ exp -^ rfcrf-- 



exp 



— ^j cxp(-2(Tfcier- — 



< 



exp I — cxp ' 



exp 



rm\' 



cxp [ —i [ T 
Sir 



Me--^ 



cxp -i rkd T ] ] ^ "^^P ^* ''"^ICI ^ 



37r 



< 



exp I — cxp ' 



e 

|2 j2 



< rk ( ^i-^ I cxp 



e 

rfcrf2 



2 cxp 



rfciei^ 



2 cxp — 



rfcieP 



<^^exp -^ +2cxp^ 



which impUes that {IV) satisfies 

v/2x 

{iv) < 



3e 



JO 

r 



r exp 



Ar'^rk 



dzdipdr 



JO 



exp dzdipdr 



< 



192 



exp 



8Tk\ ( e^h 8e 32 



eh^ J V(Tfc)2 ' Tkh ' h^J [rkf 
2tP- { Tie \^ r 



rk{h^ - 1) 



r2zE 



e^v{-y'^)dy. (13) 



10 



Also, we have the foUowing relation 



e \2 ,, ■ ■2Nfc „ / 'r'««' \ „ / ■ / , .2 37r 



2rkJ (l-^^rf) -cxp( ^lexp(-»(rfcd^- ^ 



e \ 2 



<,_J-|l-...T+exp(-lM! 



H27fc) (l + -^^T+exp 



<[^)' (l + r^ih^-im^f+cxpi- 



Tk{h^-i)\e 



which implies that 



/i /' 2 /' /i 



^0 



{l + T\h^-l)\^)''dzdipdr 
Tk{h'^ - l)r2 



cxp 



j dzdipdr 



rkh 



, v^x 



^ / (l + r^(/i^-l)^r^)'^dr 



/27r 



— / exp 



dr 






1 „ / 2Tfc(fe2_iy 

2 fh V i 



exp{~y^)dy, 
(14) 



where 



,V2^ 



'o 



T 1 

tends to zero when e < r^ and ^r- < — -. 

Taking into account the estimates obtained previously for the terms (///) 



11 



and {IV) we conclude that inside Qh we have 



dni < 



H{t) /2rfc 



(2n)i 






2h V rfc(ft2 - 1) Jo 



fv^ 



TT 

192 

27r2 



exp 



exp(-?;^)dy 



8e 32 



e/iV V(Tfc)2 rfc/i /iV (rfc) 



exp{y^)dy, 



(15) 



Taking into account the estimations obtained in (|lip and (|15l) . we conclude 
that 

\\EhA-,rk) ^ RrRhE{-,Tk)\\i^^Gn) 

Vol{Gh) 



27r 



^^(^^) i^y 



+H{rk) ('M 



exp 



h^e 



^C,{e,h,r) 



ne 



.^EK 



exp(-y^)dy 



2h V Tfc(/l2 - 1) JO 

8Tk\ / e2/j 8e 32 \ e^/i 



TT" 

192 

27r2 



exp 



exp{y'^)dy, 



2 / 2Tk 



Evaluating now the Zi— norm with respect to the time-lattice 

\\Eh.T{-,Tk) - RrRhEi- ,Tk)\\i^(G,,x(To, + oo)-,) 

= II \\Eh.r{-^Tk) - RrRhE{-,Tk)\\i^(G^)\\i^^^T^^+^)^) 

-\-OQ 

= ^ T\\EhA:^k)-RrRhE{;Tk)\\l^i^G,) 

A;— TTio + l 



(16) 



12 



VoliGh) 

27r 



+ 00 

E 

k—nio + 1 



2 e 

Trfcr 



TT 

'2/I 
192 

h 



exp 



km 
IK 



krh 



Ci{e,h,T) 



kT{h^ 



exp 



-1) 

8fcT 



hkTjh^ 



exp(-y^)dy 



e^/i 8e 32 



e^pi~y^)dy 



{kTf 



After straightforward calculations we conclude that the previous series is 
convergent and its sums, which we will denote by C2(/i, ft.,T), tends to zero 



when e < r and 



< 



1 



3.2.3 Main Result 

Using the inequalities ([9]) and (fT7|) we obtain the general estimation 



< 



To + Vol{Gh) ^K + C2(e,h,T) if r>lA(e<r8A^< gi,) 

To + Vol{Gh)^l^ + C2{e,h,T) if 0<r<lA(e<r«A^< giy) 

For the purpose of our convergence theorem we require that h < ho^ where 
ho is an arbitrary constant. Now we can formulate the following convergence 
theorem 

r2 1 

Theorem 3.1. Let —r < — ^. Then 
h'' Gtt'' 

\\E^,r - RrRkE\\,^^G,xmt) ^0 forh,T^ 0+. 

Proof. We prove that for an arbitrary S > there exists a constant h such that 

\\Eh,T - RrRhEWi^f^^^^Tg^-t-^ < S 

Also, it is possible to choose a Tq in the lattice M+, and define 

Tq^ = Tq + ar and Tq" = Tq ^ (1 ^ Oi)''' with a € [0, 1) 
such that To+ e M+ and Tq" G M+. Obviously we have 

||£;/j,^-i?r-Rh£^||,^(G^xR+) 

< \\Eh,r - RrRhE\\i^^aHyc\a.T+U + ll-^''-^ " ^^^'i^ll/i(G;,x(Ttr, + oo)0 ' 



(17) 



13 



Now, a simple estimation using ^ and ^T7\ shows that the right-hand side of 
the last inequality is bounded by 6. D 

4 Implicit Difference Equation 

In this section we will make a similar study for the implicit equation, i.e, we will 
initially obtain the fundamental solution E^ ^ of the following implicit equation 

{-AhEl,){hm,Tik + l))-t{drEl^){hm,Tk) = Sh{hm) SArk), (18) 

and then we will develop similar convergent results for this solution. 
Using the discrete Fourier transform again wc find the solution 

in analogy to Q . Finally, we get the following system of equations to calculate 

El,^{hm,Q) ^ for /im e M)^ 

{{l~iTl^h)El^A{hrn,Tk) = 5h{hrn) if /im e M^ .(19) 

((l-iTA/j)£;;;,^)(/im,T(fc + l)) = El.^{hrn,Tk) if /im e R^, rfc G Rj, fc e Z+ 

We note that it is also possible to describe the fundamental solution by 
apphcation of RhF 

ElAhm,Tk) = R^J^f^{l + tTd^)-T^A^)\{hm,Tk). 

However, this does not proves our assertion, that oi E^ ^ being a fundamental 
solution of ^. For that purpose, we need to do a similar study as in Section 
2.2. 

4.1 Existence of Fundamental Solution 

First we have the following three lemmas 

Lemma 4.1. Let fh he an arbitrary bounded function. Then the equation 

{l-iTAh)vh = fhix), 
for all X G K^ has a unique solution Vh- 

Proof. We will omit the presentation of this proof because it is very similar to 
the proof of the Lemma 1 presented in [5]. D 

The following two lemmas are proved in [5]. 

14 



Lemma 4.2. ///,, G /i(M|), then vh G h{^l)- 

Lemma 4.3. If\fh{x)\ < Xie^^il^l and \eh{x)\ < iCae"^^'^' with < ci < C2, 
then Vh{x) < K^e-'^''^-^^^''^ for all 5 > 0. 

With this three results we can present the following result about the existence 
of the fundamental solution. 

Theorem 4.4. System Iil9\) has a unique solution E*^ ^ and for arbitrary Tq < 
00 it holds El,^ e /i(M3 x [0,To]). 

Proof. The assertion follows from Lemmas 14. 1[|T^ and HT51 We remark that the 
consideration in Lemma 14.31 can be repeated as long as necessary. An estimation 
of the ?!-— norm with respect to t is possible because the number of time steps 
is bounded. D 

4.2 Convergence 

For i = we can write the difference equation in the form 

6t\ 
-i + — I El^{hrn, rk) = Shihm) 

+ -r^ [El,^{h{mi + 1), hfn2, hms, rk) + El^.^{h{mi — 1), /1TO2, hm^, rk) 

+E'^ ^{hmi, h{m2 + 1), hm^, rk) + E*^ ^{hmi, h{m2 ~ 1), /ittis, rfc) 
+£',*^(/iTOi, hm2,h{m3 + l),r/c) + E'^.^{hmi, hm2, h{m-i — l),Tfc)] , 

We have that 
6r 



-^ + 7? 



^ \ElAhm.,T)\h^ < 1 + |I ^ \El^{h7n,T)\h^ 



h 

meffif me 



which implies that 

IIKr(-,T)|L,(M3) < L 

In the same way we prove that the inequality 

IIKr(-,T(fc + i))|L,(R3) < ||i?;;^,(.,Tfc)||,^(R3), 

for each fc > 1 starting with the corresponding equations (fTOll . We obtain 

\\Eh,T\\h{GhX]0,To]^) < ll^/^,rll/i(R3x]0,To]^) 



mo 

fc=i 



15 



In relation to the continuous fundamental solution we have the result 

\\RTRhE\\h{GHX]0,To]^) 



1^- —^ 



:^ (ATrkr) 2 



VoI{Gh. 

(47rr)i T 



mo -. 



fc=i 



< 



VoUGh) -^ if T > 1 

(47r)^ 

yo/(G,,) ^i±^ if < T < 1 

(47r)2 



(20) 



Now we study 



\\ElA;Tk)-RhE{;Tk)\\i^^G,), 



para fc > ttiq- In order to estimate the right hand side of (|10p we need to study 
the following inequality, which is very similar to ^ 

\\{ThEl^)i;Tk) - iRrRhT-'E')i;Tk)\\L,(R3) 
H{Tk) (2Tk\i 



< 



(2^)^ 



2^k)' (1 



ird 



'2\-k 



exp 



Tk(f 



exp ~i 



rkS - ^ 



(V) 



+ 



exp 



Tkd' 



exp —I 



Tkd"- 



T 



exp 



—\ expl-«lrfc|C| - — 



ii(K''). 



Taking into account that {VI) = {IV) and the following relation 



{ikT (l+.rdr^^-cxpf-il^jexpf-zfrfcd^- ^ 



<(^)' (l + r2(/.,^-lf|e|Y+cxp 



rfc(/i2-l)|ep 



(21) 
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we conclude that the results obtained are equal to the conclusions obtained in 
(fT7|) . In this sense we can say that 

\\El,-RrRhE\\i,(G,^(To,+oo)^) < C2{e,KT) (22) 

T 1 

where C2(e, /i, r) is a quantity that tends to zero when e < r* and —rr < — -. 
From (HOD and dH]) we obtain 

\\Eh,T - ■Rr-R/i£'||;j(Ghx(To, + oo)) 

To + Vol{Gh) ^^ + C2{h,h,T) if r>lAe<T8A^<gi7 
To + Vol{Gh)^^^ + C2{h,h,T) if 0<r<lAe<T«A^<gi^ 

(47r) 2 " o'^ 



< 



In this conditions we can present the following theorem 
Theorem 4.5. For h ^ and t ^ we have the convergence 

\\El,{-,t) - RrR,M;mi,iG,xRt) ^ ^- 
Proof. For an arbitrary 6 we can choose h and r such that 

\\El,{;t)~RrRHE{;t)\\l,iG,.Mt) < ^- 

If we consider the same Tq of the Theorem 1 and we use Tq^ in ([20|l and T^ for 
(P^ we obtain the desired result. D 
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